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Abstract: 14-3-3a is a member of a highly conserved family of 14-3-3 proteins that has a 
double-edged sword role in human cancers. Former reports have indicated that the 14-3-3 
protein may be in an open or closed state. In this work, we found that the apo-14-3-3a is in 
an open state compared with the phosphopeptide bound 14-3 -3 a complex which is in a 
more closed state based on our 80 ns molecular dynamics (MD) simulations. The interaction 
between the two monomers of 14-3-3a in the open state is the same as that in the closed 
state. In both open and closed states, helices A to D, which are involved in dimerization, 
are stable. However, large differences are found in helices E and F. The hydrophobic contacts 
and hydrogen bonds between helices E and G in apo-14-3-3o are different from those in 
the bound 14-3-3a complex. The restrained and the mutated (Arg56 or Argl29 to alanine) 
MD simulations indicate that the conformation of four residues (Lys49, Arg56, Argl29 and 
Tyrl30) may play an important role to keep the 14-3-3a protein in an open or closed state. 
These results would be useful to evaluate the 14-3-3a protein structure-function relationship. 
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1. Introduction 

14-3-3 proteins are a family of highly conserved, acidic proteins ubiquitously expressed in all 
eukaryotic cells [1]. This family of proteins consists of seven distinct isoforms in human cells ((3, e, y, 
r|, a, t and Q as well as a variety of post-translationally modified forms [2-4]. 14-3-3 proteins act as 
adaptor proteins that control functions of its target proteins through highly regulated protein-protein 
interactions. 14-3-3 proteins have been shown to interact with several hundred proteins, some of which 
play important roles in human diseases and in the life of the cell [5]; the a isoform, for example, has 
been implicated in breast cancer [6] and is necessary for proper G2 checkpoint function [7]. Depending 
on the nature of its target protein, 14-3-3 proteins are involved in the regulation and coordination of 
many cellular processes including cell cycle progression, apoptosis, metabolism, transcriptional 
regulation of gene expression, and the DNA damage response [8,9]. 

The crystallographic structures of all seven 14-3-3 isoforms are now known [10] and have revealed 
that 14-3-3 proteins have characteristic cup-like shape functional dimers with each of the monomers 
consisting of nine antiparallel a-helices displaying a so-called amphipathic groove [10,11]. All 14-3-3 
proteins can form homodimers and heterodimers except the sigma isoform [12,13]. Each monomer can 
independently bind to a target protein in the binding groove [5] for different motifs. These motifs 
include the prototype sequence, RSxpS/TxP (mode 1), RxxxpS/TxP (mode 2) and pS/TxCOOH 
(mode 3), where x stands for any amino acid [14-16]. It seems that the 14-3-3 protein may have a rigid 
conformation because the conformation of most crystallographic structures for 14-3-3 isoforms are 
similar in both apo- and ligand-bound forms with a closed state except apo-14-3-3(3 (PDB ID: 2BQ0) 
which has been found to have one monomer in an open state while the other is in a closed state [10]. 
However, the 14-3-3 proteins can bind several hundred diverse target proteins [17,18], and a study by 
Liu et al. has suggested that 14-3-3 proteins may bind to target proteins through the transition between 
open and closed conformations [19,20]. 

The crystallographic structure (PDB ID: 1YWT) [12] reveals that the phosphorylated group of the 
binding peptide formed six hydrogen bonds with 14-3-3a protein (Figure 1). The six hydrogen bonds 
involve four hydrophilic residues (Lys49, Arg56, Argl29 and Tyrl30). The crystallographic structure 
of 14-3-3g protein (PDB ID: 1YZ5) [21] without binding peptide superimposes very well with the 
1YWT structure in the closed state. The conformation of the four residues of 1 YWT is very similar to 
1YZ5 (Figure 1). As we know, the four residues are hydrophilic and would not be stable together in 
the water environment. These reports and the crystallographic structures raise a series of interesting 
questions: (1) Is the interaction between the two monomers for the closed state different from that for 
the open state? (2) Which part of the monomer for the closed state is different from that for the open 
state? (3) Does the conformation of the four residues interacting with the phosphorylated group 
influence the state of 14-3-3a protein? To answer these questions, we utilize MD simulations to 
analyze the conformational change of 14-3-3a protein, since MD simulation is an efficient tool for 
exploring conformational changes and unfolding proteins [22-26]. 
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Here, we will discuss the equilibrium of the MD simulation, the interaction between the two monomers 
for apo-14-3-3a and bound 14-3-3a, the stability for all the helices in 14-3-3a protein and the reason 
for the conformational change to understand the differences of 14-3-3a between closed and open states. 

Figure 1. Superimposed structures of 1YWT (white) and 1YZ5 (cyan). The proteins are 
shown in cartoon representation, and the four residues and the phosphoserine are shown in 
ball and stick representations. 




2. Results and Discussion 



2.1. Equilibrium of the Molecular Dynamics Simulation 

To assess the dynamics stability of the two systems during the MD simulations, energetic and 
structural properties are monitored during the 80 ns MD simulation for both complexes. Figure 2 
shows the root mean square deviations (RMSDs) of the backbone atoms from the starting structures 
during the MD simulations. As can be seen in Figure 2, the RMSD values of each system tend to 
equilibrium after 20 ns. The RMSD values for apo-14-3-3a are larger than those for bound 14-3-3a 
during all the MD simulation time, and the averaged RMSD values for the last 60 ns MD simulation 
trajectories are 7.70 and 2.00 A for apo-14-3-3a and bound 14-3-3a, respectively. This implies that 
there is large conformational difference between apo-14-3-3a and bound 14-3-3o during the MD 
simulation process. A 60 ns MD simulation for apo-14-3-3a was carried out with different original 
speed. It has shown the similar dynamic characters as the first one. So we discuss the conformational 
change of apo-14-3-3a according to the first MD simulation. 

Figure 2. RMSDs of backbone atoms of 14-3-3 protein as a function of the MD simulation 
time. Five snapshots extracted from the first 20 ns MD simulation of apo-14-3-3a. 
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The 14-3-3 protein dimeric molecule has a characteristic cup-like shape with a central channel, and 
each monomer of 14-3-3 protein binds to the target protein by the amphipathic ligand-binding groove 
formed by helices C, E, G and I [27]. In order to understand the influence caused by the absence of the 
peptide in apo-14-3-3a, comparisons between the conformation of apo-14-3-3a and bound 14-3-3o 
were performed to obtain global information about the conformational change in the apo-14-3-3a 
structure. The last structure of bound 14-3-3o from the MD trajectories is superimposed with the 
apo-14-3-3a structure, as shown in Figure SI. For helices A, B, C and D, which are located on the 
bottom of each monomer, the structures are superimposed very well for both monomers. However, 
helices G, H and I, which are located on the top of each monomer in apo-14-3-3a, are obviously 
different from those in bound 14-3-3a. For helices E and F, which are located at the junction between 
the top and bottom of each monomer, the conformational change is not as large as the helices on the 
top. As shown in Figure SI, the conformational change can be easily detected by the distances between 
the two mass centers of helices G, H and I of the two monomers. Figure S2a plots the distances along 
the 80 ns MD simulation period for apo-14-3-3a and bound 14-3-3a. In order to evaluate the dynamic 
flexibility of helices G, H and I, the linear regressions in apo-14-3-3a and bound 14-3-3a are 
calculated for the last 60 ns and the entire MD simulation time. The distances in bound 14-3 -3a fluctuate 
around 55.74 A, as indicated by the slopes of the linear regression lines of -7.05 x 10 7 A/ps. The 
distances in apo-14-3-3a are increasing for the first 20 ns MD, and then they tend to fluctuate around 
76.53 A, as indicated by the slopes of the linear regression lines of -8.52 x 10 6 A/ps. Figure S2b and 
c show the distances between the mass center of helices A to D in the two monomers and the mass 
center of the helices G to I as a function of MD simulation time. The last structures in apo-14-3-3a and 
bound 14-3-3g are in the closed and the open states, respectively (Figure SI). The distances are less 
than 35 A in bound 14-3-3a (Figure S2b) with the closed state, as well as larger than 37 A in apo-14-3-3a 
after the first 20 ns MD simulation (Figure S2c) with the open state. The distances between the mass 
center of helices A to D in the both monomers and the mass center of helices G, H and I can be used to 
define the conformation state of 14-3-3a. 

In order to complement structural analysis, free energy calculations have been performed by using 
the MM-GBSA method. Figure 3a shows the time-series of the sums of gas-phase energies and 
solvation free energies for 14-3-3a protein without the phosphopeptide, which are computed from 
2000 snapshots extracted at 40 ps intervals from the entire MD simulations trajectories. For both cases, 
the free energies decrease at the original 20 ns MD simulation. After that, the systems tend to reach a 
stable state, the slopes of linear regression lines for the snapshots from the last 60 ns are -1.23 x 10" 4 
and -2.19 x 10~ 5 kcal/(mol.ps) for apo-14-3-3a and bound 14-3-3a, respectively. The enthalpic 
contribution is averaged over 1500 snapshots taken at 40 ps intervals from the last 60 ns of the MD 
simulations. The entropic contribution is averaged over 300 snapshots taken at 200 ps intervals from 
the last 60 ns. As shown in Table 1, the total free energy of apo-14-3-3a is lower than that of bound 
14-3-3g, this implies that the open conformation of 14-3-3a protein is more stable than the closed 
conformation. It means that peptide binding drives 14-3-3 protein transits from open to closed 
conformation. The total free energy of apo-14-3-3a protein is unfavorable for the conformational change. 
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Figure 3. (a) Sum of gas-phase and solvation free energies calculated for 2000 snapshots 
of 14-3-3g protein extracted at 40 ps intervals from the MD simulations. The linear lines 
for the last 60 ns simulation time are colored in green and blue for apo-14-3-3a and bound 
14-3-3g, respectively; (b) Root-mean-square fluctuation (RMSF) of the Ca atoms versus 
residue number for a monomer of 14-3-3a obtained from the crystallographic structure and 
the last 60 ns MD simulation. The region of helices are divided with the green vertical line 
and labeled. 




Time (ns) Residues 



Table 1. Free energies for 14-3-3o protein (kcal/mol). 



14-3-3c 


Enthalpic contribution 


Entropic contribution 


Total energy 


Apo-14-3-3o 


-5532.83 ±43.46 


2666.32 ± 6.46 


-8199.15 


Bound 14-3-3a 


-5520.95 ±43.62 


2658.52 ±6.90 


-8179.47 



The analysis for the RMSD of backbone atoms indicates that the conformational change of 
apo-14-3-3a is large. In order to investigate the stability of the nine helices composing the monomer of 
14-3-3 protein, the defined secondary structure of protein (DSSP) [28], which is commonly used to 
observe the time evolution of protein's secondary structures, is calculated. Figure S3 plots the DSSP 
along the 80 ns simulation period. There is no large differences for all the helices between apo-14-3-3o 
and bound 14-3-3a, and the helical structures are stable. However, larger differences are observed at 
the turn section, which is the junction between two helices, especially for the junction between helices G 
and H. This implies that the secondary structure of the nine helices for apo-14-3-3a and bound 14-3-3a 
is stable, the large change of apo-14-3-3a is mainly caused at the junction between the two helices. 

Figure 3b illustrates the root mean square fluctuation (RMSF) of the Ca atoms versus the residue 
number for crystallographic structure, apo-14-3-3a and bound 14-3 -3 a, respectively. As seen in 
Figure 3b, the line for crystallographic structure shows the same character as that for apo-14-3-3a or 
bound 14-3-3a. The RMSF values of the residues in the junction between the two helices are larger 
than those of the residues in the helices. Comparisons between the RMSF values for every helices in 
apo-14-3-3a and bound 14-3-3a are carried out. It is noteworthy that all the helices show the similar 
character for both systems except three helices (E, G and I). This indicates that the stabilities for the 
three helices are larger in bound 14-3-3a than in apo-14-3-3a. The crystallographic structure has 
revealed that these three helices interact with the target protein. Based on these observations, we may 
draw a conclusion that the bound 14-3-3a is more rigid than apo-14-3-3a. 
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2.2. The Interaction between the Two Monomers for Apo-l4-3-3o and Bound 14-3-3& 

As shown in Figures 2 and S2, the conformation for apo-14-3-3a continues to change and does not 
reach a stable state during the first 20 ns MD simulation. In order to investigate the interaction between 
the two monomers for apo-14-3-3a and bound 14-3-3a, the binding free energies between the two 
monomers are calculated with the MM-GBSA method. The binding free energies large fluctuation 
would indicate a large change for the interaction. Figure 4a shows the comparison of time evolution of 
the binding free energies between the two monomers in the apo-14-3-3a and the bound 14-3-3a. 
As shown in Figure 4a, the large fluctuation binding free energies for apo-14-3-3a is observed for the 
first 20 ns MD simulation; the binding free energies are stable for the both systems after the first 20 ns 
MD simulation. This is in accordance with the conformational analysis. The averaged binding free 
energies from the last 60 ns MD simulation are -25.40 and -25.21 kcal/mol for apo-14-3-3a and 
bound 14-3-3o, respectively. This implies that the dimerization affinity for apo-14-3-3a is the same as 
the affinity for bound 14-3-3a. 

Figure 4. (a) Binding free energies calculated for 2000 snapshots extracted at 40 ps intervals 
from the whole MD simulations for apo-14-3-3a and bound 14-3-3a; (b) Decomposition of 
energy on a per-residue basis into contributions from van der Waals energy CEVdw), sum of 
electrostatic energies and polar component of solvation free energy (E Q \ Q + G po i) and the 
non-polar (G non poi) component of solvation free energy for residues of apo-14-3-3a, as well 
as (c) for bound 14-3-3o; (d) Comparison between the hydrophobic interactions for per 
residue in apo-14-3-3a and bound 14-3-3a. 
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In order to further explore the detailed interaction, the binding free energies between the two 
monomers are decomposed on a per residue basis by using the MM-GBSA approach for generating a 
protein-residue interaction spectrum which is shown in Figure 4b,c. The decomposition approach is 
extremely useful to locate the residues, which contribute to the protein-protein and protein-ligand 
interactions [29-31]. The residues with the absolute free energies contribution larger than 1 kcal/mol 
for the protein-residue are labeled (Figure 4b,c). In this work, comparisons between the interaction 
spectrums of apo-14-3-3a and bound 14-3-3a are performed. The number of the favorable residues for 
dimerization is eleven and ten in bound 14-3 -3 a and apo-14-3-3o, respectively. And the number of 
the unfavorable residues for dimerization is eight and seven in bound 14-3-3a and apo-14-3-3a, 
respectively. The residues with large contribution in apo-14-3-3a are also found in bound 14-3-3a with 
large contribution. The interaction spectrum of apo-14-3-3a is very similar with that of bound 14-3-3a. 

Liu et al. have reported that the structural stability of hydrophobic cores is critical for the stable 14-3-3 
dimeric protein [19]. The contribution of the hydrophobic interaction (the sum of the van der Waals 
energy and the non-polar contribution to solvation free energy) for per residue in apo-14-3-3a are 
compared with those in bound 14-3-3a and shown in Figure 4d. The correlation coefficient r = 0.99 is 
obtained; this implies that the hydrophobic contribution for apo-14-3-3a is the same as for bound 
14-3-3g. We further determine the structural stability of the residues which have large contribution for 
the binding by calculating the RMSF. The residues unfavorable for binding are in square symbol and 
red color, and the residues favorable for binding are in circle symbol and blue color (Figure 3b). 
As shown in Figure 3b, all the RMSF values except one residue Glu75 in bound 14-3 -3a are not larger 
than those in apo-14-3-3a. This indicates that the binding pocket of apo-14-3-3a or bound 14-3-3a is 
stable, the change absent peptide in apo-14-3-3o does not cause change for the interaction between the 
two monomers. 

2.3. The Stability of the Helices A, B, C and D 

Main-chain-based clustering analysis is also performed to investigate conformational fluctuation 
and stability for helices A, B, C and D. There are six heavy atoms clusters for bound 14-3-3a, which is 
in the closed state, as well as 27 for apo-14-3-3a in the open state. The same cluster analysis is also 
performed for helices A, B, C and D; there are 5 clusters for both bound 14-3-3a and apo-14-3-3a. 
One of the five clusters is highly populated with a percentage of all the structure 49.9% and 49.8% for 
bound 14-3-3g and apo-14-3-3a, respectively. The structure which is the nearest cluster centre of the 
highly populated cluster for bound 14-3-3a is superimposed with that for apo-14-3-3a and shown in 
Figure 5a. As shown in Figure 5a, the structure of helices A, B, C and D in bound 14-3 -3 a is 
superimposed very well with that in apo-14-3-3a. The RMSD value of the backbone for the two 
nearest cluster centre of the highest populated cluster is 1.96 A. This implies that the structural changes 
of these four helices are smaller than the others in apo-14-3-3a. It is plotted for the RMSD of the 
backbone atoms of the helices A, B, C and D versus the MD simulation times shown in Figure S4a. 
Comparisons of RMSDs for the whole protein (Figure 2) and for four helices were performed. The 
RMSDs of apo-14-3-3a are higher than those of bound 14-3-3a, however, this difference is obviously 
smaller than that shown in Figure 2. 
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Figure 5. The helices of apo-14-3-3a (colored in cyan) are superimposed with the helices 
of bond 14-3-3g (colored in red) structure, (a) is for the helices A, B, C and D; (b) is for 
the helices G, H and I; (c) is for the helices E, F and G; (d) is for the helices E, F and G in 
side view with several labeled residues. 




2.4. The Stability of Helices G, H and I 

Helices G, H and I of apo-14-3-3a are obviously different from the helices in bound 14-3-3a. 
In order to investigate the conformational change and stability of these three helices, cluster analysis 
of heavy atoms was performed. There are six heavy atom clusters for the helices G, H and I for bound 
14-3-3g, as well as five clusters for apo-14-3-3a. The highest populated clusters are chosen from the 
clusters for bound 14-3-3a and apo-14-3-3a. The structure, with the smallest RMSD from the averaged 
structure in the highest populated cluster for bound 14-3-3a is superimposed with that in apo-14-3-3a 
and shown in Figure 5b. The two structures superimposed very well. This implies that the stability of 
the three helices in apo-14-3-3a is very similar to that in bound 14-3-3o. The RMSDs of the backbone 
atoms of the helices G, H and I versus the MD times were plotted in Figure S4b. 

2.5. The Conformational Change of the Helices E, F and G 

Figure S5a shows the RMSDs of the backbone atoms for the helices E, F and G from the starting 
structures during the MD simulation. The RMSDs of the helices E, F and G in apo-14-3-3a are larger 
than those in bound 14-3-3a. The structures which are the nearest to the cluster centers for the 
populated clusters of the helices E, F and G for apo-14-3-3a and for bound 14-3-3a are superimposed 
and shown in Figure 5c. It is obviously shown in Figure 5c that there are large differences for the 
relative position of the helices E and G. The distances between the mass center of the helix E and G are 
calculated during the MD simulation and shown in Figure S5b. The distances and the fluctuation of 
distance in apo-14-3-3a are larger than those in bound 14-3-3a except during first five ns MD 
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simulation. It is quite obvious that the structures of the helices E and G in bound 14-3 -3 a are different 
from those in apo-14-3-3a. 

Figure 6a shows all possible hydrophobic and hydrogen-bonding contacts between helices E and G 
in the bound 14-3-3a structure identified with LigPlot+ [32]. To further investigate the interactions 
between helices E and G, comparisons of the populations of ten hydrophobic contacts between 
apo-14-3-3a simulations and bound 14-3-3a during the last 60 ns MD simulation were performed and 
shown in Figure 6b. As can be seen in Figure 6b, the population of the ten contacts for apo-14-3-3a is 
different from that for bound 14-3-3a. Three hydrophobic contacts with populations higher than 50% 
are found in bound 14-3-3a, and no hydrophobic contacts with populations higher than 50% is found 
in apo-14-3-3a. This implies that the interaction between helices E and G for bound 14-3-3a are more 
stable than for apo-14-3-3a. 

Figure 6. (a) Two-dimensional representation for the interactions between helices E and G, 
the hydrophobic contacts and hydrogen bond are marked with red and green dash line, 
respectively; (b) The populations of ten hydrophobic contacts (labeled in Figure 9a) in the 
last 60 ns MD simulation. 




Hydrophobic contacts 

A hydrogen bond between the side chains of Asp 126 and Asnl75 is observed in the crystallographic 
structure. In order to evaluate the stability of the hydrogen bond, we investigate the hydrogen bond 
between E and G in bound 14-3-3a and apo-14-3-3a during the MD simulation. For bound 14-3-3a, 
two stable hydrogen bonds are formed between the two oxygen atoms (named OD1 and OD2) of 
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Asp 126 side chain and the nitrogen atom of Asnl75 side chain with 99.67% and 87.90% occupancy, 
respectively. However, the two stable hydrogen bonds are not found in apo-14-3-3a. In bound 14-3-3a, 
the side chain of Argl29 forms a stable hydrogen bond with the phosphorylated group. The four 
residues formed hydrogen bonds with the phosphorylated group and show stable conformation during 
the whole MD simulation (Table SI). However, a stable hydrogen bond net is found between the side 
chain of Glul82 and side chain of Argl29 in apo-14-3-3a. The minimum distances between the side 
chain of Glul82 and side chain of Argl29 are calculated for every structure during the last 60 ns MD 
simulation for apo-14-3-3a. The minimum distances are less than the largest distance 3.5 A formed, a 
hydrogen bond with 72.5%. This implies that the conformation of the four residues in the bound 
14-3-3o no longer exists in the apo-14-3-3a. 

2. 6. The Reason for the Conformational Change 

The phosphorylated group plays a key role for 14-3-3 protein to identify the target protein. To 
investigate the role of the phosphorylated group in 14-3-3a, the 60 ns MD simulation is performed for 
bound 14-3-3g without the phosphorylated group (nop-14-3-3a). The relative position of the helices A 
to D and G to I and the interactions between the two monomers are stable during the apo-14-3-3a and 
bound 14-3-3g MD simulation process. So it is reasonable to define the state of 14-3-3a by using the 
distances between the mass center of helices A to D in the two monomers and the mass center of 
helices G to I. Figure 7 plots the distances between the mass center of the helices A to D in two 
monomers and the mass center of the helices G to I tracked through MD simulations. Important 
information is found in Figure 7a, the monomers A and B are in open and closed states after the first 
20 ns MD simulation, respectively. This implies that there may be some other residues playing the 
same function as the phosphorylated group for 14-3-3a protein to bind the target proteins. 

Figure 7. Plots of distances between the mass center of the helices A to D in the two 
monomers and the mass center of the helices G to I throughout the MD simulations, (a) is 
for nop-14-3-3a; (b) is for apo-14-3-3a with four residues restrained; (c,d) are for bound 
14-3-3 with Arg56 and Argl29 mutated to alanine residues. The black and red lines are for 
the monomer A and monomer B, respectively. 
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Comparisons of the key residues in the last structure of the MD simulation for the three systems 
(bound 14-3-3g, apo-14-3-3a and nop-14-3-3a) were performed and are shown in Figure 8. It is 
obvious as seen in Figure 8a that the conformations of the four residues (Lys49, Arg56, Argl29 and 
Tyrl30) in apo-14-3-3a are different from those in bound 14-3-3a and same as those in nop-14-3-3a. 
The relative position of the peptide in nop-14-3-3a is different from that in bound 14-3-3a. The 
reasons are that the nop-14-3-3a is in an open state and the bound 14-3-3o in a closed state. This 
implies that the phosphorylated group plays an important role in the binding between the 
14-3-3g protein and the target protein. However, as shown in Figure 8b, where different conditions 
were observed, the conformations of the four residues except Lys49 in bound 14-3-3a are similar as 
those in nop-14-3-3a and different from those in apo-14-3-3a. The relative position of the peptide in 
nop-14-3-3a is similar as that in bound 14-3-3a except the C-terminal of the peptide. The C-terminal 
forms a few hydrogen bonds with the four residues. This indicates that there would be some other 
residues which play the same role as the phosphorylated residues. This result is in accordance with earlier 
studies, which show that the 14-3-3 proteins can bind their ligands in a phosphorylation-independent 
manner [27,33,34]. 

Figure 8. Superimposition of the key residues from the last structure for bound 14-3 -3a, 
apo-14-3-3a and nop-14-3-3a. The proteins are shown in cartoon representation, apo-14-3-3 
is shown in green color, bound 14-3-3 in red and nop-14-3-3a in white. The key residues 
are shown in ball and stick representation, (a) is for monomer A; (b) is for monomer B. 



As shown in Figure 8, the four residues (Lys49, Arg56, Argl29 and Tyrl30) would be in close 
proximity when 14-3-3 protein is in a closed state. Our data suggest that the state of the 14-3-3a 
protein is determined by the conformations of the key four residues. In order to test the possibility, 
another MD simulation was carried out in which the relative position of the four residues in the 
monomer A are restrained. The distances between the mass center of helices A to D in two monomers 
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and the mass center of helices G to I are plotted in Figure 7b. The conformation of monomer A is in a 
closed state during the 20 to 38 ns MD simulation times and the monomer B is in an open state. This 
implies that the conformation of the four residues may play some role for the 14-3-3a states. In order 
to evaluate the role of residues Arg56 and Argl29, two additional MD simulations for the bound 
14-3-3g with Arg56 or Argl29 mutated to alanine residues are performed. The distances between the 
mass center of helices A to D in two monomers and the mass center of helices G to I are plotted in 
Figure 7c,d. Comparisons between Figures 7c and S2b,c, reveal that the distances are larger in Figure 7 
than in Figure S2b and smaller than in Figure S2c. This implies that the residue Arg56 would play 
some function in maintaining the state of 14-3-3a protein, as well as Argl29. 

3. Materials and Methods 

3.1. System Setups 

In this study, the crystal structure of 14-3-3a determined by Wilker et al. was used as the template 
structure (PDB ID: 1YWT) [12]. Missing loops were obtained from the crystal structure of 14-3-3a 
(PDB ID: 3MHR) [5]. MD trajectories were generated for two states. One state was for bound 14-3-3a, 
the other one for apo-14-3-3a. In addition, the bound 14-3-3a systems were dephosphorylated to get 
the nop-14-3-3a systems, restrained four key residues (Lys49, Arg56, Argl29 and Tyrl30) and 
mutated the Arg56 or Argl29 to alanine, respectively. The standard AMBER force field (FF03) [35] 
was used to describe the protein parameters, the parameters of phosphoserine was designed by 
Homeye et al. [36]. The 14-3-3a proteins were solvated in a rectangular periodic box of TIP3P [37] 
water molecules with a margin distance of 12 A, and the system neutralized by adding an 
appropriate number of sodions. The apo-14-3-3a and bound 14-3-3a systems contain 74,315 and 
80,107 atoms, respectively. 

3.2. Molecular Dynamics Simulation 

All the MD simulations are carried out using the AMBER10 package [38]. Periodic boundary 
conditions and particle mesh Ewald method [39] were employed to treat long-range electrostatic 
interactions. All the covalent bonds involving hydrogen atoms were constrained by applying the 
SHAKE algorithm [40]. The integration time step for all MD simulations was set at 2 fs. The nonbonded 
cutoff was 12 A. The solvated models were first minimized with the module SANDER in constant 
volume by 1000 cycles of steepest descent minimization followed by 1000 cycles of conjugated 
gradient minimization. After energy minimization, applying harmonic restraints with force constants of 
2 kcal/(molA 2 ) to all solute atoms, canonical ensemble (NVT)-MD was carried out for 70 ps, during 
which the systems were heated from 0 to 300 K. Subsequent isothermal isobaric ensemble (NPT)-MD 
is used for 90 ps to adjust the solvent density. Finally, 80 ns isothermal isobaric ensemble (NPT)-MD 
simulation is applied to apo-14-3-3a and bound 14-3-3a simulations without any restraints. The 
nop-14-3-3a and the mutated systems are run for 40 ns. The restrained systems were run for 40 ns 
applying harmonic restraints with force constants of 2 kcal/(mol- A 2 ). The temperature was regulated at 
300 K using a Langevin thermostat and the pressure was kept at 1.0 atm using isotropic positional 
scaling. The intermediate structures were saved at every 1 ps for analysis. 
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3.3. Free Energy Calculations 

The free energies for 14-3-3a proteins were calculated by using MM-GBSA method. All water 
molecules and counterions were stripped. For every snapshot, the total energy (G tota i) was estimated 
from contribution of gas-phase energies (H gas ) and solvation free energies (G so ivation). 

total gas solvation (1) 

7/gas was further divided into a van der Waals (£Vdw) and electrostatic energies (i?eie). These energies 
were computed using the same parameters set as that used in the MD simulation. 
And the solvation free energy (G so ivation) was further divided into a polar (G po i) and a non-polar 
component (G n0 n P oi). 

^gas = ^vdW + ^ele (2) 

y^f _|_ y^r 

solvation pol nonpol (3) 

The polar contribution (G po i) was calculated by generalized Born (GB) methods implemented in 
sander. The non-polar free energy (G non poi) is determined using, 

G nonpol =YSASA+p (4) 

where SASA is the solvent- accessible surface area that is determined using the LCPO model [41]. 
The values y and (3 are the empirical constants and are set 0.005 kcal/(mol-A 2 ) and 0, respectively [42]. 

The conformational entropy contributions were estimated for 300 snapshots evenly from the last 
80 ns MD trajectories using normal-mode analysis with AMBER NMODE module [43]. 

Free energy decomposition in terms of contributions from structural subunits of both binding 
partners provides insight into the origin of binding on an atomic level. We demonstrated the decomposition 
on a per-residue basis into contributions from molecular mechanics and solvation energies but not for 
entropies. The binding interaction of each monomer-residue pair includes four terms: van der Waals 
contribution (A2? V dw) 5 electrostatic contribution (A2s e ie), polar solvation contribution (AG po i) and 
non-polar solvation contribution (AG non poi) 

AG = AE Hw + AE\ + AG , + AG , rz\ 

monmer-residue vdW ele pol nonpol yJ) 

3.4. Structural Analysis 

Ptraj [44] module of Amber Tools software is used for hydrogen bond analysis, RMSD, RMSF and 
DSSP. The formation of the hydrogen bond depended on the distance and the orientation criteria as 
follow: (1) the distance between donor (D) and acceptor (A) atoms is shorter than or equal 3.5 A; 
(2) the angle D-H...A is greater or equal to 120°. 

4. Conclusions 

MD simulations were carried out to evaluate the conformational differences between apo-14-3-3a 
and bound 14-3-3a. The apo-14-3-3a and bound 14-3-3a were in the same original closed state. 
However, apo-14-3-3a was in an open state after the first 20 ns, and bound 14-3-3a was in a closed 
state during the whole 80 ns MD simulation. Comparisons of the conformational character between the 
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open state and the closed state were performed. The relative position for the nine helices except E and 
F was stable in both systems. The interactions between the two monomers were stable for both systems. 
The conformational change for apo-14-3-3a from the closed state to the open state mainly causes 
different interactions between helices E and G. The two monomers of MD simulation for non 
phosphopeptide are in the closed and open states, respectively. This implies the 14-3-3a protein may bind 
to not only phosphopeptide but also non phosphopeptide. The restrained four residues (Lys49, Arg56, 
Argl29 and Tyrl30) MD simulation indicates that the conformation of these four residues play an 
important role for 14-3-3a protein in an open or closed state. The four could bind to not only the 
phosphorylated residues but also other molecules with stronger polarity. These results would be helpful to 
understand the structure- function relationship of 14-3 -3 a protein, and to support drug design of 
14-3-3g protein. 
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